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I. 


INTRODUCTION 


The  Influence  of  pressure  on  the  burning  rates  of  propellants  is  of 
great  interest  to  ballisticians.  No  adequate  theoretical  model  now 
exists  to  explain  the  experimental  facts,  such  as  the  familiar  power  law. 
One  approach  to  this  problem  is  to  investigate  the  pressure  dependence 
of  the  activation  energies  of  chemical  reactions.  The  rationale  is  made 
clear  by  Figure  I,  which  shows  the  potential -energy  curves  (approximated 
by  the  harmonic  oscillator)  along  some  configuration  coordinate,  Q,  for 
the  ground  (lower  curve)  and  an  excited  (upper  curve)  electronic  state 
of  a molecule.  The  activation  energy  is  that  which  the  curves  cross  on 
the  left.  Clearly,  pressure  effects  on  the  shapes  and  positions  of 
these  curves  can  change  the  activation  energy.  An  attempt  to  derive 
such  a pressure  dependence  of  activation  energy  from  a published  model 
of  pressure  effects  on  electronic  spectra  revealed  serious  deficiencies 
in  the  sK>del.  This  report  covers  a study  of  these  deficiencies  and  a 
possible  remedy. 

1 2 

Drickamer  e^  ^ ’ developed  a simple  and  attractive  model  for  the 
effects  of  high  pressure  on  molecular,  electronic  spectra.  Okamoto 

3 4 5 

^ ' ' proposed  and  carried  out  experimental  tests  of  the  model. 

Lin^  has  provided  a more  rigorous  theoretical  foundation  for  the  simple 
model,  generalizing  it  in  the  process.  I shall  demonstrate  that  the 
experimental  tests  are  inconsistent  with  the  Drickamer  model,  but  not 
with  the  Lin  theory. 

Part  II  summarizes  the  Drickamer  model  and  the  Okamoto  tests. 

Part  III  exhibits  fallacies  in  the  authors'  data  analysis  and  shows  that 
the  data,  in  fact,  contradict  the  model.  Part  IV  summarizes  the  Lin 
theory  and  casts  the  equations  in  a form  amenable  to  experimental  test. 
Part  V reinterprets  Okamoto* s data  in  the  context  of  Lin's  theory.  The 
results  are  discussed  in  Part  VI,  and  concluding  remarks  follow  in 
Part  VII. 


n Drickamer,  H.G.,  Frank,  C.W. , and  Slichter,  C.P.,  Proc.  Nat.  Acad. 
Sci.  (USA)  933  (1972). 

2.  Drickamer,  H.G.,  and  Frank,  C.W.,  "Electronic  Transitions  and  the 

High  Pressure  Chemistry  and  Physics  of  Solids,"  London:  Chapman 

and  Hall  (1973). 

3.  Okamoto,  B.Y.,  Drotning,  W.E.,  and  Drickamer,  H.G.,  Proc.  Nat.  Acad. 
Sci.  (USA)  71,  2671  (1974.  Note  that  R is  misprinted  R^  in  the 

last  term  o?~their  Eq.  (1).  The  correct  equations  appears  in  ref.  4. 

4.  Okamoto,  B.Y. , and  Drickamer,  H.G.,  J.  Chem.  Phys.  2870  (1974). 

5.  Okamoto,  B.Y.,  Thesis,  Univ.  of  Illinois,  Urbana-Champaign  (1975). 

6.  Lin,  S.H. , J.  Chem.  Phys.  59,  4458  (1973). 
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II.  THE  DRICKAMER  MODEL 

The  basic  model  assumes  that  pressure  couples  to  harmonic  oscilla- 
tions along  a single  configuration  coordinate,  viz.,  volume.  The 
potential  functions  for  the  ground  and  excited  electronic  states  are 
then  given,  respectively,  by 


V ■ ♦ Qp 


(2.1) 


Eq  ♦ W^(Q-q)^  ♦ Qp. 


(2.2) 


where  u is  the  force  constant,  p is  pressure,  Q is  the  configuration 
coordinate,  and  Eq  and  q are  the  energy  and  coordinate  separations, 
respectively,  between  the  potential  minima.  Projection  of  one  potential 
minimum  upon  the  other  potential  function  in  accordance  with  the 
Franck-Condon  principle  yields  expressions  for  the  optical  absorption 
and  emission  energies,  respectively: 


* E + JsRu^q^  + Rqp  + >s(R-l)P^/u^ 


(2.3) 


Eq  - ♦ qp/R  ♦ ’s(R-l)P  /(R  u ), 


(2.4) 


where  R Similar  projection  of  the  terminal  points  of  oscillation 

of  the  lowest  vibrational  level  gives  the  Gaussian  half-widths. 


N |Ruq  * (R-l)p/ai| 


(2.5) 


c - N|<u  q/R**  ♦ (R-l)p/(r^/2a.)f  , 


(2.6) 


where  N > (8kT  log  2)^.  (References  1 and  2 also  generalize  the  model 
to  multiple  configuration  coordinates  and  normal  modes,  but  only  the 
basic  form  was  used  by  Okamoto.) 

3 

Okamoto  ^ al  elaborated  upon  Eqs.  (2.3)  - (2.6)  and  proposed  a 
procedure  for  esHmatlng  the  unknown  parameters  from  spectral  data. 
They  expressed  the  peak  shifts  and  broadening  by 


«|  ■ Rqp  ♦ (R-l)p^/2«^  ♦ %Rq^(o)^-uJ)  , 


(2.7) 
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(2.8) 

(2.9) 

(2.10) 


, » qp/R  ♦ (R-1)p^/2rV  - >jq^(o.^-«J), 

■ n(|  Rqw  ♦ (R-l)p/u»|  - I Rq«Q  |), 

«Ce  ■ n(I  R"/*^**  * (R-l)p/R^^^uil  - Iqu^l  ) , 

2 2 
where  u is  assuned  to  be  pressure  dependent,  and  Ug  is  the  value  at 

atoospheric  pressure.  The  pressure  dependence  of  R is  assumed  to  be 
relatively  insignificant. 

4 

Okamoto  and  Drickamer  experimentally  examined  the  effects  of 
pressures  up  to  8 GPa  on  the  absorption  and  emission  spectra  of  anthra- 
cene, phenanthrene,  and  tetracene,  both  as  neat  crystals  and  in  various 
solvents.  They  found  that  the  data  could  be  fitted  by 

«^a  ■ *iP  * (2.11) 

and 

- d^P.  (2.12) 

along  with  the  corresponding  expressions  with  e substituted  for  a. 

Okamoto' s thesis^  provides  details  of  the  data  analysis  and  lists  the 
experimental  parameters.  These  have  been  converted  to  SI  units  and 
displayed  in  Table  I.  In  order  to  evaluate  the  parameters  of  Eqs.  (2.7) 

- (2.10),  Okamoto  made  the  further  assumption, 

(1)^  ■ u)q(1  ♦ bjp  + b2P^),  (2.13) 

where  b.  and  b,  are  empirical  parameters  derived  from  the  compressibility 

* ^ 7 2 

data  of  Vaidya  and  Kennedy  . Thence  he  estimated  R,  q,  and  and 
concluded  that  the  model  was  internally  consistent. 


III.  FALLACIES  IN  THE  OKAMOTO- DRICKAMER  RESULTS 

Attempts  to  apply  these  results  in  other  equations  of  the  model 

revealed  serious  discrepancies.  Examination  of  Okamoto 's  thesis^  shows 
that  upon  substitution  of  Eq.  (2.13)  into  Eqs.  (2.7)  - (2.10),  certain 
approximations  were  made.  For  example,  Eq.  (2.7)  became 

■ RqP  ♦ (R-1)P^/2«5  ♦ «|Rq^»J(bjP  ♦ b2P^)  . (3.1) 


T.  Vaidya,  S.N.,  and  Kennedy,  G.C.,  J.  Chem.  Phys.  987  (1971). 
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In  essence  he  assuned  that  in  the  second  term  on  the  right  hand  side, 

(bj  ♦ b2P)p  « 1.  (3.2) 

Inasmuch  as  insertion  of  the  smallest  values  of  b^  and  b2  used  gives  a 

value  of  2.3  at  only  2 GPa,  the  assumption  is  hardly  justified.  Further- 
more, the  approximation  is  unnecessary.  The  algebra  for  the  full 
expression  is  tedious  but  not  difficult. 

Another  approach,  however,  is  more  revealing.  It  will  be 
convenient  first  to  circumvent  the  annoying  absolute  values  in 
Eqs.  (2.9)  and  (2.10).  We  note  that 

.2.2  r ^ - o , . ..  .2  , o,.  . 

fie  = e - (e  ) = (fie)  +2  ce-(c)  = (fie)  + 2e  (fie), 

where  e**  is  the  half-width  at  atmospheric  pressure.  Eqs.  (2.9)  and 
(2.10)  are  now  replaced  by 

fief/N^  » rV(<*>^  -“«)  ♦ 2R(R-l)qp  + (R-l)^p^/oi^  (3.4) 

and 

fie^/N^  - q^(o)^  - u^)/R  ♦ 2(R-l)qp/R^  + (R-1) . (3.5) 

For  Eq.  (2.12)  we  write 

fieJ/N^  - 2E^p  ♦ D^p^,  (3.6) 

2 2 o 

where  B • a/w  and  E ■ d e , with  a similar  expression  for  emission. 

Aft  AAA 

(This  formulation  has  the  additional  advantage  of  utilizing  more  of  the 
experimental  data.)  Values  for  D , E , D , and  E are  given  in  Table  I. 

A A C C 

We  obtain  from  Eqs.  (2.7),  (2.11),  (3.4),  and  (3.6),  after  simpli- 
fication, 

Rq^(u)^)^  ♦ 2{Rq-aj)p  - 2a2P^  - Rq^wJ  u)^  + (R-l)P^  = 0 (3.7) 

and 

2222  22572  22  (3.8) 

R“q‘‘(u^)'  ♦ 2R(R-l)qp  - 2E^p  - D^p‘‘  - R q (u^  u ♦ (R-1)  p * 0 

2 

These  are  two  quadratic  equations  in  u having  the  same  roots  for  any 
given  value  of  p.  From  the  relations  between  roots  and  coefficients 

R(R-1)Vp^-  R^(R-l)qV.  (3.9) 
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whence  R = I,  Making  that  substitution  in  Eqs.  (3.7)  and  (3,8),  we 
obtain 

q^(ta)^  - 0)^)  2(q  - a,)p  - 2^2?^  ~ 0 (3.10) 

and 


q^u)^  - 0)^)  - 2E^p  - D^p^  = 0.  (3.11) 

2 

Having  made  no  assumption  as  to  the  functional  dependence  of  u upon  p, 
we  have  obtained  expressions  having  the  form  of  Eq.  (2.13)!  Similar 
analysis  for  emission  yields 

q^((D^  - u)^)  + 2(e^  - p)p  + 2©2P^  = 0 (3.12) 


and 

q^(u)2  - u^)  - 2EgP  - Oy  = 0. 
We  see  immediately  that 


(3.13) 


(3.14) 

(3.15) 


and 


*1  - ®1  * ^a  * ^e  • 

These  relations  are  not  supported  by  the  data  in  Table  I. 

We  must  conclude,  therefore,  that  the  experimental  data  of  Okamoto 

4 

and  Drickamer  do  not  support  the  simple  model  of  Reference  2. 


IV.  THE  THEORY  OF  LIN 


Lin  considered  a solute  molecule  dissolved  in  rigid  solvent 
molecules.  He  wrote  a pressure-dependent  Hamiltonian  for  the  electronic 
motion  and  obtained  the  potential  function  for  nuclear  motion  and  the 
electronic  wave  function  by  a first-order  perturbation.  His  potential 
functions  for  ground  and  excited  states  are 

V * Ejj  ♦ E(*iU)jQj  ♦ CjQjP)  (4.1) 

and 


11 


(4.2) 


W-^  (Q.-q.)^  c'.Q.p 
j Kr 


E'  > E 
o 

J 

where  j references  the  different  normal  modes  and  coordinates,  and  the 
Cj  are  coupling  constants.  The  energy  difference  between  the  potential 

minima  (referred  to  by  Drickamer  as  the  thermal  energy)  is  then 


E c'jqjP  + >i(cj/a)j 


/2/  72,  2 

c'  /(/  )p 


(4.3) 


where  E..  (0)  » E' ^ - E . 
th  00 

These  expressions  are  the  same  as  those  given  by  Drickamer  for  the 
general  case  with  one  exception,  pointed  out  by  Lin'  Drickamer  implicitly 

assumed  d.  * c..  Lin  also  noted  another  difference  from  Drickamer's 
J J 

treatment,  which  does  not  show  up  in  these  equations.  First-order  per- 
turbation theory  makes  the  normal  frequencies  (w^  etc.)  independent  of 

pressure.  In  what  follows  I shall  follow  Lin  on  these  points. 


Whereas  Drickamer  appealed  to  simple  geometric  arguments  to  derive 
his  expressions  for  the  spectral  shifts,  Lin  carried  through  a full 
treatment  of  the  band  shape  function.  For  the  absorption  spectra  he 
obtains  the  following  expressions,  analogous  to  Eqs.  (2.3)  and  (2.5), 

£ = E..  + I ? - hE(a)'./2(Dj)(u.-(i)' )(3u.-oj'. )coth(ho)./2kT)  (4.4) 

a tn  J .jjjjj  Jj  J 

and 

■ 4 log  2 
a 

♦ (b^/u»j)(u)^-o)')^coth^(ha)j/2kT)  j , (4.5) 

where 

- qj  ♦ (Cj/«oj  - (4.6) 

and  E^j^  is  given  by  Eq.  (4.3).  The  corresponding  expressions  for  emission 

are  obtained  by  interchanging  Wj  and  wherever  they  appear  explicitly  in 

Eqs.  (4.4)  and  (4.5)  and  by  changing  the  signs  of  the  second  and  third 
terms  on  the  right-hand  side  of  Eq.  (4.4). 


E j}  rh(a)'^/o)J6  ^coth(hio./2kT) 
j J L J J J 


Some  simplification  is  necessary  before  we  can  make  practical  use  of 
these  equations.  Let  us,  then,  make  the  following  two  assumptions: 
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1.  The  force  constants  for  excited  and  ground  states  stand  in 
the  same  ratio  in  all  normal  modes  and  configuration  coordinates,  i.e., 

2 2 

Rj  * 01^  /(iit  can  be  replaced  by  R. 

2.  The  coupling  constants  for  excited  and  ground  states  stand  in 

the  same  ratio  in  all  normal  modes  and  configuration  coordinates,  i.e., 
rj  “ replaced  by  r. 

We  can  now  factor  R and  r out  of  the  summations  and  define  the  new 
parameters , 


and 


q « I c.q. 

j 


1/u  ^ * E c?/(i)^  . 

j ^ ^ 


(4.7) 


(4.8) 


Remembering  that,  for  small  x,  coth  x 1/x  and  making  the  appropriate 
substitutions,  we  obtain 


- 2-2  2 
6S^  * ’s(R-2r+l)p  /(i)  =>  ajP  + 82?  , 

6ef/N^  a 2R(R-r)qp  + (R-r)  V/«  ^ - 2E  p + D p^, 

fit  A 


(4.9) 

(4.10) 

6B^  » (Rr+r-R)qp/R  + >sr (2R-Rr-r)p^/R^  w ^ = e^p  + e2P^  , (4.11) 


and 


«cW  » 2(R-r)qpV  + (R-r)V/R^  w ^ = 2E  o * D . (4.12) 

V.  REINTERPRETATION  OF  THE  OKAMOTO'S  DATA  IN  CONTEXT  OF  LIN'S  THEORY 
We  have  immediately  from  Eqs.  (4.9)  - (4.12),  for  absorption. 


and 


Rq  » aj, 

R - 2r  1 a 2a2  u 
R(R  - r)q  ■ E^, 


(R  - r)^  - 


(5.1a) 

(5.2a) 

(5.3a) 

(5.4a) 
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► 


and 

for  emission. 

(Rr  ♦ r - R)q  » e^R, 

(5.1e) 

r(2R  - Rr  - r)  ■ e2R^  u ^ > 

(5.2fc) 

and 

(R  - r)q  - E^R^, 

(5.3e) 

2 3-2 

(R  - r)  - D R"*  u)  ^ . 
© 

(5.4e) 

Of  the  aany  possible  combinations  of  these  eight  equations  to  estimate 

the  four  parameters,  R,  r,  q,  and  u from  the  data  of  Table  I,  five  seemed 
most  straightforward  and  meaningful. 

Method  1 uses  only  absorption  data,  Eqs.  (5.1a)  - (5.4a).  Method  2 
uses  only  emission  data,  Eqs.  (5.1e)  - (5.4e).  Method  3 uses  only  the 
shifts  in  the  peaks,  Eqs.  (5.1a),  (5.1e),  (5.2a),  and  (5.2e).  Note  that 
either  pair  of  equations,  (5.3)  or  (5.4),  yields  R directly.  Method  4 
uses  Eqs.  (5.4a)  and  (5.4e),  complemented  by  (5.1a)  and  (5.1e).  Method 
5 uses  Eqs.  (5.3a)  and  (5.3e),  complemented  by  (5.1a)  and  (5.2a).  The 
explicit  expressions  for  the  parameters  are  displayed  in  the  Appendix. 

4 5 

Neither  Okamoto  and  Drickamer  nor  Okamoto  reveals  the  experimental 

g 

temperature,  required  to  compute  N.  However,  Mitchell  ^ al  have 
recently  reported  similar  experiments  in  the  same  laboratory  carried 
out  at  297K.  I have  assumed  this  temperature  in  the  calculations. 

Table  II  lists  the  results  of  the  calculations. 


j VI.  DISCUSSION 

The  results  can  be  evaluated  only  in  terms  of  internal  consistency. 
In  general  the  best  consistency  appears  in  hexane  solution,  the  poorest 
in  the  neat  crystals.  The  most  consistent  parameter  is  q,  the  pressure 
shift  of  the  configuration  coordinate;  the  least  consistent  is  u>,  the 
nonaal  frequency.  The  poor  results  for  u can  be  attributed  to  small 
differences  between  large  numbers.  The  variability  in  R and  r suggests 
that  my  simplifying  assumptions  may  be  too  gross. 

! 4 

t The  experimenters  questioned  all  the  emission  datA  for  crystalline 

anthracene  and  tetracene,  and  some  of  those  for  crystalline  phenanthrene 
because  of  possible  exciaer  formation.  It  is  surprising,  therefore, 
that  Method  1,  which  uses  only  absorption  data,  fails  for  crystalline 


S.  Mitchell,  D.J.,  Schuster,  G.I.,  and  Drickamer,  H.G.,  J.  Am. 
Chem.  Soc.  99,  114S  (1977). 
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TABLE  II.  Configuration-Coordinate  Parameters 


System 

R 

r 

q 

(cm^  mol”^) 

(lO^Pa.j'^'^^) 

Anthracene  crystal 

Method  1 

-1.09 

-.164 

6.5 

2.3 

Method  2 

.55 

.69 

-14.6 

.78 

Method  3 

.65 

.89 

-10.9 

b 

Method  4 

.96 

1.40 

- 7.4 

1.08 

Method  5 

.99 

1.92 

-14.2 

b 

Phenanthrene  crystal 

Method  1 

.25 

.59 

-12.2 

1.82 

Method  2 

.65 

.78 

- 4.6 

1.15 

Method  3 

.58 

.69 

- 5.3 

3.0 

Method  4 

.90 

1.09 

- 3.4 

1.10 

Method  S 

.90 

1.23 

- 3.6 

b 

Tetracene  crystal® 

Method  1 

- .060 

.41 

170.3 

1.52 

Method  2 

.76 

.89 

-13.3 

b 

Anthracene  in  PIMA^ 

Method  1 

.153 

.35 

-20.0 

1.67 

Method  2 

.74 

.84 

- 3.8 

1.23 

Method  3 

.87 

.94 

- 3.5 

b 

Method  4 

.99 

1.08 

- 3.1 

.78 

Method  5 

.98 

1.17 

- 3.3 

b 

Phenanthrene  in  PfWA*^ 

Method  1 

.75 

.90 

- 1.21 

1.67 

Method  2 

-1.30 

- .22 

.75 

b 

Method  3 

.83 

.92 

- 1.10 

.90 

Method  4 

.88 

.98 

- 1.04 

2.1 

Method  5 

.71 

.86 

- 1.56 

1.57 

Anthracene  in  hexane 

Method  1 

.82 

.85 

-10.8 

.40 

Method  2 

.76 

.78 

-11.6 

.54 

Method  3 

.94 

.96 

- 9.5 

.133 

Method  4 

1.08 

1.11 

- 8.2 

.39 

Method  5 

.83 

.85 

-13.6 

.40 

Phenathrene  in  hexane 

Method  1 

.69 

.75 

- 3.9 

.94 

Method  2 

.76 

.80 

- 3.6 

.92 

Method  3 

.84 

.89 

- 3.2 

.48 

Method  4 

.95 

1.02 

- 2.8 

1.08 

Method  5 

.85 

.91 

- 4.3 

.36 

*Data  not  available  for  Methods  2,  4 and  S. 
^laaginary  value. 

^olyaethyl  aethacrylate  fila. 


IS 


anthracene  and  tetracene  (negative  R)  and  appears  to  be  out  of  line 
for  crystalline  phenanthrene  and  for  anthracene  in  PhMA.  The  poor 
performance  of  Method  1 is  the  most  puzzling  feature  of  the  results. 

The  theory  of  Lin^  promises  to  provide  a useful  approach  to  the 
study  of  pressure  effects  on  the  electronic  structure  of  complex 
molecules.  It  appears  already  to  be  a useful  qualitative  tool,  but 
refinements  will  be  needed  to  obtain  a formulation  combining  adequate 
simplicity  and  accuracy  for  practical,  quantitative  applications.  ^ 
Much  more  experimental  data  similar  to  that  of  Okamoto  and  Drickamer 
will  be  required  for  a good  evaluation. 


VII.  CONCLUSIONS 

4 

A reexamination  of  the  experimental  tests  by  Okamoto  and  Drickamer  ’ 

1 2 

of  the  simple  Drickamer  model  * of  high-pressure  effects  on  electronic 
spectra  revealed  errors  in  the  data  analysis.  Properly  interpreted,  the 

^ • 

tests  fail  to  verify  the  model.  A theory  of  Lin  was  found  to  be  more 
compatible  with  the  data,  but  it  will  require  further  refinement  before 
it  is  fully  adequate  for  practical  applications. 

It  should  now  be  possible  to  take  a preliminary  look  at  pressure 
effects  on  activation  energy  and  perhaps  reaction  rates,  in  a quali- 
tative way.  Such  work  is  now  in  progress. 
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APPENDIX 


This  Appendix  gives  the  explicit  expressions  for  estimating  the 
configuration-coordinate  parameters  by  the  five  methods  of  Part  V. 

Method  1 


(R  - r)  - E^aj 

(A.l) 

- R-r)^/D^, 

(A.2) 

r ■ 1 - 2^2^^  * (R-r) . 

(A. 3) 

R - (R-r)  ♦ r. 

(A. 4) 

q « a^/R. 

(A. 5) 

Method  2 

(l-o^R^  ♦ R^  ♦ Oj(2-Oj)R  ♦ Oj  = 0, 

(A.  6) 

r - R(R+Oj)/(R(R+l)  ♦ o^)  , 

(A. 7) 

q - E^R^/(R-r), 

(A. 8) 

- (R-r)^/D  R^, 

(A. 9) 

where  Oj  • e^/E^,  O2  “ Zej/D^. 

Method  3 

R^  - (l-P2)/(p2  - 2pjP2  +pj). 

(A. 10) 

r - R(pjR+l)/(R+l), 

(A. 11) 

q ” a^/R, 

(A. 12) 

■ (R  - 2r  ♦ l)/2a2. 

(A. 13) 

where  p^  ■ (>2  “ ®2^*2‘ 

Method  4 

r’  ■ d,/d,. 

(A. 14) 

q ■ ®j/R, 

(A. 15) 

r • R(ej^q)/q(R«l), 

(A. 16) 

■ (R-r)^ZD^. 

(A. 17) 
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Method  S 


«’  • E./E,. 

(A. 18) 

(R-r)  ■ E,/*,, 

(A. 19) 

r ■ R - (R-r). 

(A. 20) 

q - E^/R(R-r), 

(A. 21) 

■ (R  - 2r  ♦ l)/2a2  • 

(A. 22) 
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